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We search for the production of doubly charmed baryons in e + e~ annihilations at or near a center- 
of-mass energy of 10.58 GeV, in a data sample with an integrated luminosity of 232 fb -1 recorded 
with the BABAR detector at the PEP-II storage ring at the Stanford Linear Accelerator Center. We 
search for S^ c baryons in the final states A£ K~tv + and Ecir + , and S^ c + baryons in the final states 
A+ K~ 7r + 7r + and S^tt + -k + . We find no evidence for the production of doubly charmed baryons. 

PACS numbers: 13.85.Rm, 14.20.Lq 



I. INTRODUCTION AND OVERVIEW 

The lowest-mass doubly charmed baryons are pre- 
dicted to be the members of an isospin doublet (5+. = ccd 
and = ecu [1]) with J p = \ and L = 0. There 

are many theoretical predictions for the 3~£ c and 
masses and lifetimes [2-11]. The predicted masses lie in 
a range of approximately 3.5 to 3.8GeV/c 2 [3-8]. The 
mass difference between the S+ and the is pre- 

dicted to be on the order of lMeV/c 2 [9]. The H+ 
and Sj c + lifetimes are expected to be between about 
0.1 and 0.2ps, and 0.5 and 1.5ps, respectively [10, 11]. 
Theoretical estimates for branching fractions relevant to 
this paper are B(S+ -»■ A+A"~7r+) = 0.03, B(S+ -> 
S° c tt+) = 0.02, B(E++ -> A+K-TT+TT+) = 0.05, and 
B(S++ - H c °7r+7r+) = 0.05 [12]. 

Several predictions have been made for the production 
cross sections of doubly charmed baryons in e + e~ anni- 
hilations [13-15]; the predictions range from 1 to 250 fb 
for an e + e~ center-of-mass (CM) energy near 10.58 GeV, 
and translate into 0(1O 2 -1O 4 ) doubly charmed baryons 
produced in the BABAR data set of 232 fb -1 analyzed 
here. Measured cross sections for double-cc production 
in Belle [16] and BABAR [17] are an order of magnitude 
larger than non-relativistic QCD predictions. Calcula- 
tions for cccc and cccc cross sections are very similar; 
therefore, the predicted cross sections for doubly charmed 
baryons may also have been underestimated. 

The SELEX collaboration, which uses the Fermilab 
600-GeV/c charged hyperon beam, has published evi- 
dence for the baryon in the A+K~ir + and pD + K~ 
decay modes with a mass of (3518.7 ± 1.7) MeV/c 2 [18, 
19]. The £"+ + baryon, detected in the decay mode 
A+ K~ir + Tr + , with a mass of 3460MeV/c 2 , was reported 
by SELEX at ICHEP 2002 [20]. The 5+-S++ mass dif- 
ference of 60McV/c 2 is not consistent with theoretical 
expectations. SELEX sets an upper limit (at 90% con- 
fidence level) of 33 fs on the lifetime of the E+ c baryon, 
in conflict with theoretical predictions. The photopro- 
duction experiment FOCUS does not observe any E cc 
states [21] although they observe 19, 500 A+ baryons, 
compared to 1, 650 for SELEX. 



*Also at Laboratoirc de Physique Corpusculairc, Clermont- 
Ferrand, France 

tAlso with Universita di Perugia, Dipartimcnto di Fisica, Perugia, 
Italy 

t Also with Universita dclla Basilicata, Potenza, Italy 



In this paper, we describe a search for the production 
of S cc baryons in a data sample corresponding to an inte- 
grated luminosity of 232 fb -1 recorded with the BABAR 
detector at the PEP-II asymmetric-energy e + e~ storage 
ring at the Stanford Linear Accelerator Center. Events 
containing A+ — ► pK~ir + candidates are searched for the 
presence of £"+. — > A+K~n + and H+ + — > A+K~n + n + 
candidates. Events containing A — ► pir~ candidates are 
searched for the presence of — > S®ir + and S+ c + — > 
E®ir + ir + candidates where E~ir + and E~ — > Air~ . 

The BABAR detector is described in detail else- 
where [22]. The tracking of charged particles is pro- 
vided by a five-layer double-sided silicon vertex tracker 
(SVT) and a 40-layer drift chamber (DCH). Discrimina- 
tion among charged pions, kaons, and protons relies on 
ionization energy loss (dE/dx) in the DCH and SVT, 
and on Cherenkov photons detected in a ring-imaging 
detector (DIRC). A CsI(Tl) crystal calorimeter is used 
to identify electrons and photons. These four detector 
subsystems are mounted inside a 1.5-T solcnoidal super- 
conducting magnet. The instrumented flux return for the 
solcnoidal magnet provides muon identification. 

For event simulations, we use the Monte Carlo (MC) 
generators JETSET74 [23] and EVTGEN [24] with a full 
detector simulation based on GEANT4 [25] . These simu- 
lations are used to estimate the reconstruction efficiencies 
of the searches. For each of the four S cc decay channels 
used in our searches, we produce approximately 100,000 
simulated e + e~ — » cc events in which at least one of 
the primary charm quarks hadronizes into a E cc . The 
distribution of momentum in the CM frame (p*) for sim- 
ulated S cc peaks at about 2.5GeV/c, with 80% above 
2.0GcV/c and 62% above 2.3GcV/c. The 5+ and S++ 
baryons are simulated with the SELEX masses of 3520 
and 3460MeV/c 2 , respectively. The E+, S++ , and A+ 
decays are generated according to phase space. 

We search for E cc production as an excess of candi- 
dates in the distribution of the difference in the mea- 
sured masses of the S cc and the candidate daughter 
baryon. Some mass uncertainties cancel in this mass dif- 
ference, improving the mass resolution and thereby the 
signal-to-background ratio. We use the following nota- 
tion: AM (A - B) = M(A) - M(B), where A is the 
parent and B is the daughter baryon. M(X) refers to 
the measured invariant mass of the X candidate. 

Selection criteria are chosen to maximize e/s/B, where 
e is the simulated reconstruction efficiency and B is the 
number of candidates in data in upper and lower side- 
bands of the mass-difference regions in which we search 
for S cc signals. During this process the search regions 
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were hidden to minimize potential experimenter bias. 

Charm hadrons carry a significant fraction of the ini- 
tial energy of the charm quark, whereas random combi- 
nations of charged particles in an event form lower-energy 
candidates. To take advantage of this difference, we se- 
lect S cc candidates for which the p* of the S cc is above 
a minimum value. For S cc decay modes containing a 
A+, the optimal requirement is p* > 2.3GcV/c. Because 
the background levels for events containing a S c candi- 
date are lower, we apply the less stringent requirement 
p* > 2.0GeV/c. To facilitate comparisons with theoret- 
ical predictions, we repeat the searches with no require- 
ment on p* . 

We conduct searches for E cc near the masses of the 
states observed by SELEX and over wider ranges that 
include many of the theoretically predicted masses. We 
use MC techniques to account for the width of the search 
region in the statistical interpretation of the results. 



II. SEARCH FOR DECAYS TO A+K--k + (tt+) 

In the searches for E*. — > A+K~tt + and S+ + — * 
A+ K~tt + tt + , we reconstruct the A+ baryon in its decay 
to pK~ 7r + . Pion, kaon and proton candidates are identi- 
fied using the SVT, DCH and DIRC. The x 2 probability 
for the A+ daughter particles and for the E cc daughter 
particles to each come from a common vertex is required 
to be above 1%. The number of reconstructed At signal 
events is approximately 600, 000. 

The distribution of the mass difference AM (E cc — 
At) is shown in Fig. 1 for candidates with M(A+) 
between 2281 and 2291MeV/c 2 (±0.8cr), and also for 
M(A+) sidebands (2256 < M(A+) < 2281MeV/c 2 and 
2291 < M(A+) < 2316MeV/c 2 ). To search for a 
signal in data and to estimate the efficiency, we per- 
form two-dimensional fits to M{At) and AM(E CC — 
A+). The range of M{A+) used in all fits is 2256 to 
2316 McV/c 2 . We search for E cc states with masses be- 
tween 3390 and 3600MeV/c 2 (AM(S CC - A+) between 
1100 and 1310 MeV/c 2 ). The mass-difference sidebands 
in data are between 890 and 1100 MeV/c 2 , and 1310 and 
1520MeV/c 2 . 

Approximately half of all background E cc candidates 
are due to true A+ particles combined with random 
pion and kaon candidates from the rest of the event. 
This background is fit with a Gaussian shape in M(A+) 
and a linear shape in AM(E CC — Another signif- 

icant background contribution is from false Aj candi- 
dates. This source of background is fit with the product 
of a linear function in M(At) and a linear function in 
AM(S CC -A+). 

MC simulations show that n cc signals peak in three dif- 
ferent ways in the M{A+) versus AM(E CC — At) plane. 
In most cases, the E cc is reconstructed correctly and the 
measured values of both M(A+) and AM(E CC — yl+) lie 
close to the generated values; such candidates are fit with 
the product of two Gaussian distributions, one in each 




1.2 1.3 1.4 1.5 

AM(S -At) (GeV/c 2 ) 



FIG. 1: Distributions of the mass difference AM(S CC — At) 
for (a,b) Etc and (c,d) Ett~ candidates with (a,c) no p* re- 
quirement and (b,d) p* > 2.3 GeV/c. Data points with 
error bars correspond to candidates near the At mass: 
2281 < M(At) < 2291MeV/c 2 . Shaded histograms corre- 
spond to candidates in M(vl+) sidebands (2256 < M(At) < 
2281MeV/c 2 and 2291 < M(A+) < 2316MeV/c 2 ), scaled to 
represent the expected amount of non-At background in the 
data projections. 



variable. The MC signal resolution for AM(^ CC — A+) 
is 3.5MeV/c 2 and 3.0McV/c 2 for S+ and »++, respec- 
tively. When S cc candidates are reconstructed from the 
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correct tracks but the kaon and/or pion from the 
decay is swapped with the kaon and/or pion from the 
E cc decay, the reconstruction has the correct M(S CC ) 
but an incorrect M(A+). These events are fit in both 
MC simulations and data with a Gaussian function in 
AM(S CC - A+) + M(A+) = M(S CC ) and are included 
as part of the signal. When the yl+ is correctly recon- 
structed but is combined with an incorrect pion and/or 
kaon to form the S cc , the reconstruction has the correct 
M(A+) but an incorrect AM(S CC - Af). Such events are 
not distinguishable from Aj combinatoric background. 

Each shape parameter describing the signal is con- 
strained in the fit to lie within a range determined from 
the Monte Carlo simulation, allowing for possible inac- 
curacies in the simulation. The integral of the signal 
function is allowed to be negative. Efficiencies for the 
reconstruction of E cc baryons decaying to Aj K~it + and 
A+ K~ir + Tr + are calculated from the signal yields from 
fits to the MC simulated samples. These efficiencies are 
listed in Table I. The systematic uncertainties are due 
to inaccuracies in the simulation of tracking reconstruc- 
tion (0.8% per track, added linearly) and particle iden- 
tification (1.0% per kaon, 1.0% per pion, and 4.0% per 
proton). When setting upper limits on production cross 
sections, additional systematic uncertainties arise due to 
the uncertainties on the integrated luminosity (1.0%) and 
a(e+e- -> A+X)B(A+ -> pK~x + ) (4.7%). 

TABLE I: Efficiencies determined from e + e~ — > S CC X sim- 
ulations. With the p* criterion applied, the efficiency is cal- 
culated for E cc baryons generated with p* above 2.3 GeV/c 
for the modes and 2.0 GeV/c for the S c modes. The first 
error is statistical; the second is systematic. 



p* Criterion Particle A+ Mode Eff. (%) Mode Eff. (%) 

Yes 4.2 ±0.1 ±0.2 6.1 ± 0.1 ± 0.6 

Yes S+ 10.4 ±0.1 ±0.5 9.3 ± 0.1 ± 0.7 

No S+ + 3.6 ±0.1 ±0.2 5.9 ±0.1 ±0.5 

No S+ 9.7 ±0.1 ±0.5 9.0 ±0.1 ±0.7 



We conduct searches for a signal within 10-MeV/c 2 - 
wide regions around the £"+ and masses reported 

by SELEX, and within the 210-MeV/c 2 -wide region de- 
scribed earlier. The wide search region is divided into 
21 sequential 10-MeV/c 2 search sub-regions. For each 
sub-region, we perform a two-dimensional fit over a 100- 
MeV/c 2 -wide range in mass difference centered on the 
sub-region, constraining the mean of the Gaussian signal 
function to lie within that sub-region. 

The significance of any potential signal is determined 
through the use of parametrized MC simulations. Sam- 
ples of pairs of variables (M(A£), AM(S CC — A+)) are 
generated according to the background shapes measured 
in data, with no signal contribution. The distributions of 
M(yl+) versus AM(S CC — A+) from these simulations are 
then searched in the same manner as in data. A signif- 
icance measure N/on, where N is the fitted number of 
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FIG. 2: Distributions of the mass difference AM(E CC -E°) for 
(a,b) S^ c and (c,d) S^ + candidates with (a,c) no p* require- 
ment and (b,d) p* > 2.0 GeV/c. Data points with error bars 
correspond to S cc candidates reconstructed using E° candi- 
dates near the S° mass, 2457 < M(S°) < 2485MeV/c 2 ; the 
shaded histograms correspond to M(E°) sidebands (2451 < 
M{~°) < 2457MeV/c 2 and 2487 < M{5°) < 2501MeV/c 2 ) 
scaled to represent the expected amount of non-£'° back- 
ground in the M(3°) signal region. 



signal candidates and on is the uncertainty on this num- 
ber, is determined for each fit. fn order to statistically 
combine the results of the 21 fits into one search, only 
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the largest of the 21 significance measures is used. The 
significance measure from data is compared to the distri- 
bution of significance measures from the MC simulations 
that represent those data. This comparison gives the 
probability of measuring this particular value of TV/er/v 
or higher in data under the hypothesis that no S cc are 
produced. 

None of the A+ decay mode searches finds evidence 
for E cc . The most statistically significant signal is for 
a S+ baryon with AM(S+ - A+) between 1250MeV/c 2 
and 1260MeV/c 2 , when candidates are required to have 
p* > 2.3GeV/c. With a significance measure of N/<jn = 
66/24, we find that there is an 8% probability that back- 
ground alone could produce this signal. This corresponds 
to a significance of 1.4 a, which does not constitute evi- 
dence for the 3~£ c baryon. 

Using efficiencies (e) listed in Table I and integrated lu- 
minosity (L) of (232 ± 2) fb _1 , we extract values for the 
upper limit on the production cross section times branch- 
ing fraction(s) (iS) directly from negative-log-likelihood 
functions. A conversion factor T — Le and its uncer- 
tainty ajr are incorporated in a Gaussian extension to 
the likelihood function (£) so that all systematic uncer- 
tainties are included in the results. C takes the form 

<.r-tf N 

c = e -( N -Sf-n b ) e -TSr- p( . . ^ L n ^ d) ^ 

i 

where TV is the total number of fitted events; Sf = n s 
and rib are the fitted number of signal and background 
events, respectively; / is the fitted conversion factor from 
S to n s ; a are shape parameters; and P is the probability 
function for the data point xt. The value of iS for which 
— ln£ is 1.35 units above the minimum value for which 
S is positive is interpreted as the 95%-confidence-level 
upper limit. These limits are listed in Table II. 

To facilitate comparison with the production rate of 
and to take advantage of the cancellation of the — > 
pK~ir + branching fraction, we also normalize the upper 
limits to <r(e + e~ — > X)B(A^ — > pK~ir + ), measured 
with 22fb _1 of data collected at yfs ~ 10.54 GeV; these 
upper limits arc also listed in Table II. The p* criterion 
that is applied to the S cc candidates is also applied to 
the A+ candidates in the normalization mode. 

III. SEARCH FOR DECAYS TO S°7r+(7r+) 

In the search for E£ c — > 5®tt + and S+ + — > 5®tt + it + 
decays, the is detected in the decay chain 5® — > 
S~n + , S~ — > Air~ , A — > pn~ . We search for E cc states 
with masses between 3370 and 3770MeV/c 2 {AM(S CC - 
E°) between 900 and 1300 MeV/c 2 ). The mass- difference 
sidebands in data are 800 < AM(S CC -S°) < 900 MeV/c 2 
and 1300 < AM(S CC - S°) < 1400MeV/c 2 . 

For A and 5~ candidates, we require a minimum 
signed three-dimensional flight distance of +2.0 cm and 
+0.5 cm, respectively, where the flight distance is the pro- 
jection of the vector from the primary vertex to the decay 



point, onto the momemtum vector of the candidate. A 
candidates are required to be within ±3.6MeV/c 2 (±3er) 
of the world average mass [26]. E~ candidates are re- 
quired to be within ±5.4MeV/c 2 (±3tr) of the world av- 
erage mass difference AM(S" — A), and E% candidates 
are required to be within ±14 MeV/c 2 (±2cr) of the world 
average mass difference AM(H C ° — S~) [26]. For all can- 
didate baryons, we require the vertex fit to have a \ 2 
probability greater than 0.01%. The number of recon- 
structed H° signal events is approximately 11, 700. Fig- 
ure 2 shows the distributions of mass difference for all 
S cc candidates that satisfy these criteria, with no p* re- 
quirement and with p* > 2.0 GeV/c. The reconstruction 
efficiencies are given in Table I. 

Systematic uncertainties arise mainly from possible in- 
accuracies in the simulation of track reconstruction and 
particle identification (5% for and 6% for ver- 
tex quality (6%), and mass and mass-difference resolu- 
tions (1%); the values in parentheses are the relative un- 
certainties in these efficiencies. Other sources include 
uncertainties in the total luminosity (1.0%) and in the 
branching fractions for A — > pn~ (0.8%) and S~ — > Air~ 
(0.03%). 

To search for a signal in the 400-MeV/c 2 -wide search 
region, we fit the mass-difference distribution with two 
Gaussian functions, with common means and fixed 
widths, to represent the signal, and a first-order poly- 
nomial for the background. The values of the Gaussian 
widths are determined from the MC simulation; the root- 
mean-squared deviation for AM(S , + — 5®) is 5.5 MeV/c 2 
and for AM(H++ - E°) it is 4.2 MeV/c 2 . We conduct 
50 fits with the mean of the Gaussian signal function 
constrained to lie in 50 10-MeV/c 2 ranges, each of which 
overlaps neighboring ranges by 2 MeV/c 2 . Using a MC 
approach, we calculate the upper limit on the number of 
signal events using the statistically most significant of the 
50 fits. To do this, we generate N signal events according 
to the Gaussian signal function and background events 
according to a first-order polynomial, where the num- 
ber of background events is determined from the mass- 
difference sidebands. We fit the resulting MC distribu- 
tion as described above for data, and record the number 
of signal events S for the statistically most significant fit. 
We repeat this process 10,000 times, varying TV by the 
fractional systematic uncertainty on efficiency. We then 
find the value F for which only 5% of the trials have 
S < F. We repeat the above process starting with dif- 
ferent values of TV to find the value of N for which F 
is the number of signal events found in the most signif- 
icant fit in data. This value of TV is the 95% CL up- 
per limit on the number of events, shown in Table II for 
both and , with and without p* requirements. 
We also present in Table II the limits obtained when we 
explicitly search for the states observed by SELEX. For 
comparison, the measured rate for the singly charmed S c 
baryon in BABAR is a{e+er -> S°X)B(»° ~> ~-tt+) = 
(388 + 39 + 41) fb [27]. 
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TABLE II: The 95%-confidence-level upper limits on measured rates for the production of S cc baryons with and without a 
p* requirement of 2.3 GeV/c for A + modes and 2.0 GeV/c for the S c modes. The columns labeled A r+ ' + - ) give the upper 
limits on the number of signal 3£ c ' baryons. a +< - + ' > denotes the production cross section a(e + e~ — > Etc X); & in the 
denominator indicates that the cross section has been normalized to a(e + e~ — > A+ X)B(A + — ► pK~n + ). The factor B in a 
column heading signifies that the values in the column correspond to a cross section times the branching fractions B(Et} + ^ — * 
A+K--K + (iY + ))B(At -> pK~-K + ) for decay modes with At and B(E+ (+) -> S°n + {■k+))B{S° -> S~Tr + ) for decay modes with 
5°. For each wide mass range, the upper limit corresponds to the maximum upper limit over the range. 







Upper Limits for 


S cc ^ A?K-tv + (tv + ) 


Upper Limits for 








N+ 


a + B (a + /a)B 


N++ a++B (a++/a)B 


N+ 


a+B 


JV++ 


a++B 


Wide Mass Range 


328 


14.5 fb 13.3 x 10" 4 


199 23.9 fb 22.0 x 10~ 4 


58 


4.3 fb 


58 


6.6 fb 


Wide Mass Range, p* Req. 


106 


4.4 fb 5.6 x 10~ 4 


54 5.5 fb 6.9 x 10~ 4 


41 


3.0 fb 


28 


3.1 fb 


SELEX Mass 


169 


7.5 fb 6.9 x 10~ 4 


91 10.9 fb 10.0 x 10~ 4 


26 


2.0 fb 


49 


5.6 fb 


SELEX Mass, p* Req. 


53 


2.2 fb 2.7 x 10" 4 


31 3.2 fb 4.0 x 10 -4 


18 


1.3 fb 


31 


3.4 fb 



IV. SUMMARY 

In conclusion, we have searched for doubly charmed 
baryons in e + e~ annihilations at or near a center-of-mass 
energy of 10.58 GeV. We do not observe any significant 
signals for the baryon in the decay modes A + K~w + 
and £®ir + , or for the S+ c + baryon in the decay modes 
A+K~7r + 7r+ and S°tt + tt+. 
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